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The Weyl semimetal (WSM) is a newly proposed quantum state of matter. It has Weyl nodes in 
bulk excitations and Fermi arcs surface states. We study the effects of disorder and localization in 
WSMs and find three exotic phase transitions. (I) Two Weyl nodes near the Brillouin zone bound¬ 
ary can be annihilated pairwise by disorder scattering, resulting in the opening of a topologically 
nontrivial gap and a transition from a WSM to a three-dimensional (3D) quantum anomalous Hall 
state. (II) When the two Weyl nodes are well separated in momentum space, the emergent bulk 
extended states can give rise to a direct transition from a WSM to a 3D diffusive anomalous Hall 
metal. (Ill) Two Weyl nodes can emerge near the zone center when an insulating gap closes with 
increasing disorder, enabling a direct transition from a normal band insulator to a WSM. We deter¬ 
mine the phase diagram by numerically computing the localization length and the Hall conductivity, 
and propose that the exotic phase transitions can be realized on a photonic lattice. 

PACS numbers: 72.15.Rn, 73.20.Fz, 73.21.-b, 73.43.-f 


Topological quantum states of matter have emerged as 
an important and growing field in condensed matter and 
materials physics recently 0,i- The Weyl semimetal 
(WSM) is a newly proposed quantum state of the kind 
that breaks time-reversal symmetry or inversion symme¬ 
try 00]. A WSM exhibits a set of paired zero-energy 
Weyl nodes (linearly touching points of conduction and 
valence bands) in its bulk spectrum and Fermi arcs ex¬ 
citations localized on the surface. A number of candi¬ 
date materials have been predicted to be WSMs, includ¬ 
ing pyrochlore iridates and magnetic topological insula¬ 
tor multilayers Recently, following the theoretical 

prediction [7|, angle resolved photoemission experiments 
confirmed that TaAs is a WSM by the observation of 
the Fermi arcs surface states 0. Both the Weyl nodes 
and the Fermi arcs have been observed in NbAs using 
a combination of soft X-ray and ultraviolet photoemis¬ 
sion experiments 0. Furthermore, the Weyl points have 
been predicted and subsequently observed remarkably in 
gyroid photonic crystals [101 ] . 

In this Letter, we study both numerically and analyt¬ 
ically the stability of the gapless Weyl nodes and Fermi 
arcs against random potential scattering and the novel 
disorder-induced metal-insulator transitions in WSM sys¬ 
tems. Previous studies have concentrated on the prop¬ 
erties of a single Weyl node, assuming that the disor¬ 
der potential is smooth enough to avoid scattering be¬ 
tween different nodes [Tl0T3 ]. Indeed, a system with a 
single Weyl node is not subject to Anderson localiza¬ 
tion even for strong disorder [l5[. However, the theorem 
of Nielsen and Ninomiya states that gapless Wev l nodes 
with opposite chirality must appear in pairs [161 ]. Thus, 
it is essential to study the localization properties of a 
pair of Weyl nodes since they can be annihilated pair¬ 


wise when approaching each other in momentum space 
or by strong intervalley scattering 0,0. To this end, 
we study a model system of quasi-3D WSMs formed by 
stacking 2D quantum anomalous Hall (QAH) layers along 
the z-direction Q which can support a single-pair of Weyl 
nodes in the first Brillouin zone (BZ), and uncover a rich 
set of novel disorder and localization effects. Consider 
the case when the Weyl nodes are located close to each 
other, which can happen either near the zone boundary 
or the zone center. We find that in the former case, the 
Weyl points are unstable against disorder, which causes 
the pair to merge at the zone boundary and annihilate 
each other, whereby opening up a topologically nontriv¬ 
ial gap in the bulk. This gives rise to a transition from 
the WSM to the 3D QAH state with fully quantized Hall 
conductivity. In the latter case, we find that the criti¬ 
cal state with two overlapping Weyl nodes at the zone 
center is stable against disorder as they would repel each 
other and lift the degeneracy in favor of a WSM. More¬ 
over, if one starts with a normal insulator with a band 
gap at the zone center, increasing disorder would lead to 
the closing of the band gap and the emergences of two 
repelling Weyl nodes near the zone center, i.e. a nor¬ 
mal band insulator to WSM transition. It is thus rather 
remarkable that tuning the disorder strength can switch 
the system between a WSM and a normal band insulator, 
as well as between a 3D QAH and a WSM, in analogy to 
the “phase switching” between a normal band insulator 
and a topological insulator by disorder 0-22]. A pair of 
Weyl nodes well separated in momentum space are found 
to be stable against weak disorder. However, we find that 
strong enough disorder induces bulk extended states and 
gives rise to a direct transition from the WSM to a 3D 
diffusive anomalous Hall metal 
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This disorder 
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induced transition is unconventional since it takes place 
between two metallic states and is only enabled by the 
topological character of the WSM. We study the various 
multiple disorder-driven phase transitions and obtain the 
phase diagram by evaluating the localization length and 
the Hall conductivity 
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We also performed cal¬ 
culations using the self-consistent Born approximation 
(SCBA) in order to gain important analytical insights 
into the stability of the Weyl nodes at weak disorder. 
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FIG. 2: (Color online). The phase diagram on the W — m z 
plane. The symbols guided by the solid lines are obtained 
from the localization length. The two solid triangles corre¬ 
spond to the phase boundaries in the clean limit. The blue 
dashed lines are the phase boundaries determined using the 
SCBA. 


FIG. 1: (Color online). The phase diagram of the Weyl 
semimetal Hamiltonian Ho on the t z /mo-m z /mo plane. The 
open and filled black circles correspond to the parameter val¬ 
ues studied for the disordered phase diagram shown in Fig[2] 
and the filled black circles to the localization length and Hall 
conductivity shown in Figs [3] and [4] 


We begin with the simple 2x2 WSM Hamiltonian writ¬ 
ten in momentum space [4] 

Hq = ( m z — t z cos k z )(j z + mo(2 — cos k x — cos k y )a z 
+t x a x sin k x + t y a y sin k v , (1) 


with the model parameters t x , y , z , m z , mo and Pauli ma¬ 
trices cr x , y ,z- The lattice wave vectors fc x>y . z are defined 
in the first BZ of a L x x L y x L z cubic lattice with the 
lattice constant a = 1. H 0 can be regarded as describing 
a layered (labeled by k z ) Dirac system coupled along the 
^-direction by the inter-layer coupling t z . Diagonalizing 
Hq gives the energy-momentum relation 


E k = ±yj A 2 + (t x sin k x ) 2 + (t y sin k y ) 2 (2) 


where A- = m z — t z cos k z -I-mo(2 — cos k x — cos k v ). Note 
that Hq breaks the time reversal symmetry when the first 
two terms (proportional to a z ) in Eq.Jl]) are nonzero, i.e. 
when A 2 ^ 0. Adding the random on-site potential to 
Hq, we obtain the full Hamiltonian with disorder 


H = H 0 + 


bi(r) 0 \ 

0 V 2 (r) J 


(3) 


where V) 2 (r) are uniformly distributed within 
\—W/2,W/2] with W representing the disorder strength. 


The phase diagram of the clean WSM Hamiltonian Hq 
is shown in FiglU When t z = 0, the layers decouple. 
There are three critical points on the phase diagram de¬ 
noted by the red solid squares with massless Dirac nodes 
located at the 2D BZ center (0,0) for m z /mo = 0 and 
at the zone boundaries (0,7r) and ( 7 r, 0) for m~/mo = —2 
and (-7T, 7r) for m./mo = —4. They separate three gapped 
phases: the normal insulator (NI) and two QAH phases 
with 2D Hall conductances a^ y = — e 2 //i (left) and e 2 //i 
(right) respectively [30]. Turning on the inter-layer cou¬ 
pling t z naturally leads to two 3D QAH phases with a 
bulk topological gap and fully quantized Hall conduc¬ 
tivity a xy = ±e 2 //i @. More importantly, Fig Jl] shows 
that the WSM phase emerges between the 3D QAH 
states and the NIs and is bordered by a pair of phase 
boundaries emanating from each of the three 2D critical 
points described by, from left to right, m z = ±f z — 4mo, 
±t z — 2mo, and ±t z . The WSM phase supports pairs 
of gapless Weyl nodes split by t z along the fc z -direction 
and are located at [tt, 7r, ± arccos(m 2 + 4mo0)/f z ]; 
[0,7r, ± arccos(m 2 + 2mo)/t z ] and [7r, 0, ± arccos(m z + 
2mo)/t z ]; and [0, 0, ± arccosm 2 /t z ] between the corre¬ 
sponding pairs of phase boundaries respectively. 

In the rest of the paper, we will, for the sake of simplic¬ 
ity, focus on the parameter region —0.714 < m 2 /mo < 
1.429 in the clean phase diagram (Fig. [T]) by setting 
t x , y ,z/‘mo = 0.476 and investigate the disorder induced 
phase behaviors with varying disorder strength W. The 
WSM phase in this regime (—1 < m 2 /t z < 1) has a 
single pair of Weyl nodes located at (0,0, ±fc°), k° z = 
arccos (m z /t z ). They move from the zone center (0,0,0) 
to the zone boundary (0, 0, it) as the WSM traverses the 
NI and the 3D QAH phase boundaries. 
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Our main results are summarized in the phase diagram 
in the presence of disorder, which is shown in Figj2] on 
the W—m, plane. There are four distinct phases and five 
possible multiple phase transitions with increasing disor¬ 
der strength W: (i) QAH -metal-NI, (ii) WSM-QAH 
metal-NI, (iii) WSM-metal-NI, (iv) NI-WSM-metal- 
NI, and (v) NI-metal-NI. The phase boundaries indi¬ 
cated by the symbols and guided by the solid lines are 
obtained by numerical computations of the localization 
length and the Hall conductivity. Note that the WSM- 
QAH and WSM-NI phase boundaries at weak disorder 
agree remarkably well with the ones obtained analytically 
using the SCBA described in the supplemental section. 
The findings of such weak disorder induced transitions 
between electronic states of different topological charac¬ 
ters broaden the concept of topological Anderson insula¬ 
tors U 221 and support the generality of the rich physics 
behind the interplay of randomness and topology. 



FIG. 3: (Color online), (a)-(f) Normalized localization length 
A = \(L)/L versus disorder strength W at different mass 
m z . The curves correspond to different sample widths L = 
8,10,12,14. Other parameters are mo = 2.1 and t Xl y lZ = 1. 

To calculate the localization length, we consider a 3D 
long bar sample of length L z and widths L x = L y = 
L with periodic boundary conditions in the x and y- 
directions. The localization length A (L) of the sample is 
calculated using the transfer matrix method |24l - |26| . In 
general, the normalized localization length A = X(L)/L 
increases with L in a metallic phase, decreases with L in 
an insulating phase, and is independent of L at the crit¬ 
ical point of the phase transition. The scale-dependent 
behaviors of A versus the disorder strength W, obtained 
at different values of m z and shown in Fig. (3[ reveal a 
sequence of disorder induced Anderson transitions in our 
model. In Figj3])a), the decrease of A with increasing L 
at small W places the system at m z /mo = —0.714 in 
the 3D QAH state at weak disorder on the phase dia¬ 
gram (Fig0. Increasing W encounters a critical point 
at W ~ 6.0 where dA/dL = 0, beyond which a diffusive 
metallic phase emerges as indicated by dA/dL > 0. The 


latter has a nonzero Hall conductivity (see below) and is 
thus identified with an anomalous Hall (AH) metal. Fur¬ 
ther increasing W causes another metal-insulator transi¬ 
tion in Fig(31)a) as the system becomes an Anderson insu¬ 
lator at strong disorder. Such disorder-induced multiple 
phase transitions (QAH - AH metal-NI) are analogous to 
those found in multilayer quantum Hall structures and 
3D Z 2 topological insulators 32-341]. 

Repeating this procedure at different values of m z /mo 
maps out the disordered phase diagram shown in Fig [2] 
Note that the phase boundary between the QAH and 
WSM is slanted toward the WSM side, which indi¬ 
cates that a clean WSM in proximity to the QAH state 
(i.e. when the Weyl nodes are close to = n) would 
be unstable against disorder W and would undergo a 
sequence of WSM-QAH metal-NI transitions with in¬ 
creasing W. This unexpected behavior can be deduced 
directly from the scale-dependence of the localization 
length in Figs[3]]b)-(c) obtained at m z /mo = —0.452 and 
—0.381. To understand this result, we performed SCBA 
valid for weak disorder analytically in the supplemen¬ 
tal section [3f|, which amounts to determine the band 
renormalization induced by disorder [Hill. Interest¬ 
ingly, the most important disorder-induced renormaliza¬ 
tion in SCBA is for the topological mass m z , which en¬ 
hances the band inversion and causes the Weyl nodes to 
“attract” each other to the zone boundary. When the 
they meet at the zone boundary, the Weyl nodes annihi¬ 
late pairwise and open up a nontrivial bulk gap, leading 
to the transition from the WSM phase to the 3D QAH 
phase. This phase transition can also be characterized by 
the disappears of the Fermi arcs in favor of fully quan¬ 
tized Hall conductivity to be discussed later. One can 
apply this analysis to the phase boundary between the 
WSM and the NI in Figj2] and find another surprising 
result. The WSM in proximity to the NI with the pair of 
Weyl nodes close to the zone center is stable against weak 
disorder since the effect of renormalizing m z by disorder 
will cause the Weyl nodes to repel and move apart. This 
explains why the phase boundary is slanted toward the 
NI side in Fig|2]and thus an unexpected transition from 
a NI to a WSM with increasing disorder. Physically, as 
the insulating band gap closes, a pair gapless Weyl nodes 
is nucleated at the zone center. The scale-dependence of 
A shown in Fig|3[f) at m z /mo = 0.571 confirms this with 
the switching from dA/dL < 0 to dA/dL > 0 at a criti¬ 
cal W ~ 3.3. Thus, a pair of Weyl nodes with opposite 
chirality are stable against intervalley scattering at weak 
disorder, provided that they are sufficiently away from 
the zone boundary. Remarkably, for strong enough disor¬ 
der, we find a direct transition from the WSM to the AH 
metal that enables the WSM-metal-NI and NI-WSM- 
metal-NI transitions shown in the phase diagram. Such a 
phase transition comes from the emer gen t bulk extended 
states induced by strong disorder [lj, |23| and is thus an 
unusual transition between two metallic states as can be 
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FIG. 4: (Color online), (a)-(f) Ensemble averaged Hall con¬ 
ductivity a xy as a function of the disorder strength W on 
30 x 30 x 8 samples. The parameters are from m z = —1.5 to 
rriz = 1.2 in one to one correspondence to Fig[3]]a)-(f). The 
red dash-dotted lines are guide to the eye. The error bars are 
magnified ten times to show conductance fluctuations. 


seen from the positive scale-dependence dA/dL > 0 on 
both sides of the critical point in Figs[3]]d)-(e), and on the 
two sides of the second scale-invariant point in Fig|3]T). 
Indeed, the WSM and the 3D diffusive AH metal can 
be distinguished by the nature of the Hall conductivity 
quantization. 

To study the Hall conductivity and the Fermi arcs in 
WSMs, we view the Hamiltonian Hq as a set of 2D mas¬ 
sive Dirac Hamiltonian with fixed k z as in Burkov et 
al. [BJ. The 2D quantized Hall conductance in the x — y 
olane at a fixed k z is given by a xy {k z ) = 0(fc° — \k z \)e 2 /h 
1,0, 30, 31]. The total Hall conductivity of the system 


k z eBZ a ayy /where k z takes on quantized 


iS (7 X y — 

values in the first BZ of L z layers under the periodic 
boundary condition. Thus a xy is proportional to 2 kP z , 
the distance between the two Weyl nodes. It increases 
from zero to e 2 /h as the pair of Weyl nodes moves from 
k z = 0 at the WSM-NI boundary to k z = ±tt at the 
WSM-3D QAH boundary. As a result, the Hall conduc¬ 
tivity can characterize the positions of the Weyl nodes as 
well as the density of Fermi arcs surface states. 

In the presence of disorder, we calculate the ensemble 
averaged Hall conductivity a xy of quasi-3D samples with 
L x x L y x L z = 30 x 30 x 8 using the noncommutative 
Kubo formula under periodic boundary conditions in all 
three directions [271 - 1291 . 1.3(| . The results are shown in 
FigsQIa)-(f) at the corresponding parameters of the lo¬ 
calization lengths in Figs[3]]a)-(f). Overall, a xy increases 
with the disorder strength W initially and then decreases 
with W in the strong disorder limit. FigHKa) shows that 
the 3D QAH phase sustains a fully quantized conduc¬ 
tivity a xy = e 2 //i until the system enters the diffusive 


AH metal phase at W ~ 6.0 where bulk extended states 
emerge and cr xy looses quantization, decreases with W, 
and exhibits conductance fluctuations, in complete agree¬ 
ment with the phase diagram obtained from the localiza¬ 
tion length shown in Fig[3]]a). FigsQJb)-(c) start in the 
WSM phase close to the QAH phase boundary and thus 
show the “fractional” quantized a xy = 7e 2 /8 h. Upon in¬ 
creasing W, the pairs of Weyl nodes approaches the zone 
boundary and annihilate each other at a critical W, while 
the Fermi arcs turn into a set of chiral edge states and 
the Hall conductivity acquires the maximum quantized 
value e 2 //i of the QAH state [37]. Further increasing W 
leads to the sequence of QAH-metal-NI transitions. In¬ 
terestingly, Figs[U]d)-(e) show that when the Weyl nodes 
are sufficiently away from the zone boundary, increas¬ 
ing disorder leads to the discretized growth of the dis¬ 
tance between the Weyl points and thus the increase in 
the fractions of the quantized Hall conductivity. How¬ 
ever, the Weyl nodes do not reach the zone boundary 
before being destroyed by disorder amid a direct transi¬ 
tion from WSM to 3D diffusive AH metal as the Fermi 
arcs scatter through the emergent bulk extended states 


14 1231 ] . The Hall conductivity is no longer quantized. It 


decreases with disorder and vanishes at the transition to 
the Anderson insulator at strong disorder. Thus, a Weyl 
semimetal can be distinguished from a diffusive anoma¬ 
lous Hall metal by the quantized fractions of the Hall 
conductivity. Finally, Fig|4jf) shows that as the band 
gap closes with increasing disorder, a NI with <r xy = 0 
can make a transition to a WSM at a critical W ~ 3.3, by 
nucleating a pair of Weyl nodes at the zone center such 
that a xy increases from 0 to e 2 /8 h and then to 3e 2 /8 h 
as the distance between the nodes and the density of the 
Fermi arcs states increase with disorder. At W ~ 8 and 
before the Weyl nodes can reach the zone boundary, the 
WSM-AH metal transition takes place. We stress that 
all the phase behaviors obtained from the Hall conductiv¬ 
ity are in quantitative agreement with those determined 
from the localization length in Fig [3] and displayed in 
the phase diagram in FigEJ including the location of the 
phase boundary. This further demonstrate the consis¬ 
tency and reliability of the obtained results. 


In summary, we studied the effects of disorder and An¬ 
derson transition in the simplest WSMs with two Weyl 
nodes and obtained an unexpectedly rich phase diagram. 
We find that weak disorder has important effects when 
the Weyl nodes are close in momentum space. The 
pair-annihilation near the zone boundary and the pair- 
nucleation at the zone center leads to WSM to QAH as 
well as NI to WSM transitions. Moderately strong disor¬ 
der, on the other hand, produces a WSM to diffusive AH 
metal transition. These results have important implica¬ 
tions for real WSM materials discovered recently [Ml 
which have multiple pairs of Weyl nodes. Furthermore, 
the WSM phase has been observed recently in gyroid 
photonic crystals where disorder can be introduced and 
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controlled by a speckled beam E30. In photonic crys¬ 
tals, the bulk gap can be detected by bulk transmission 
measurements while the highly directional surface trans¬ 
mission spectroscopy can be used to detect the chiral sur¬ 
face state through its unidirectional group velocity Q. 
We therefore propose that the exotic phase transitions 
found here can be readily tested on optical lattices such 
as in gyroid photonic crystals. 
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I. SELF-CONSISTENCE BORN APPROXIMATION (SCBA) 


In this section, we provide some details of the SCBA calculation. We follow the discussions of Groth et. al.— who 
performed SCBA calculations for the weak disorder induced transition from a normal band insulator to a topological 
insulator in systems with time reversal symmetry. We start with the clean part of the Hamiltonian for WSMs, i.e. 
Eq.(l) in the main text, 

Hq = ( m z — t z cos k z )a z + mo(2 — cos k x — cos k y )a z + t x a x sin k x + t y a y sin k y , (1) 


where the model parameters t XyV ^ z , m z , m q have the same meanings as in the main text. Adding the random on-site 
potential to Hq gives the full Hamiltonian in the presence of disorder 


H = H 0 + 


Vi(r) 0 \ 

0 V 2 (r) J 


where V\^{p) are uniformly distributed within [— W/2,W/2\ with W representing the disorder strength. 


( 2 ) 



Figure S 1: (Color online). The phase diagram on the W — m z plane reproduced partly from Fig.2 in the main text. The 
scattered symbols guided by the solid lines are obtained from finite size scaling of the localization length. The blue dashed 
lines are the phase boundaries obtained using the SCBA. The solid dots and squares connected by the vertical dashed lines are 
the cuts where the SCBA renormalized band dispersions are shown in Fig.S [2] They are labeled by the roman numerals (I) to 
(VI) in one to one correspondence to the six vertical panels in Fig.S O 

Within the single-particle picture, the effects of disorder can be accounted for in terms of the self-energy defined as 

- - - = ( - - - ) 

E f -H q -E(E f ) ' E f — H 


( 3 ) 
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where (...) denotes the average over disorder configurations. It is convenient to expand the 2x2 self-energy matrix 
in the basis of Pauli matrices 

= So<70 + fCTy + S z <r z , (4) 

which allows us to identify the renormalized topological mass and the chemical potential from the self-energy correc¬ 
tions, 


m z = m z + lim Re E z , 

k—> 0 

Ep = Ep + lim ReTjQ. 

/ c—>-0 


(5) 
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Figure S 2: (Color online). Panels (I)-(VI) show the evolution of the SCBA dispersion with increasing disorder strength W 
(moving up vertically) and increasing mass m z (moving to the right horizontally). The dispersion is calculated from a 3D long 
bar sample of length L z = oo along the ^-direction and widths L x = L y = 14a under periodic boundary conditions in the x 
and y directions. The SCBA renormalized mass fh z is obtained by solving Eq. [5] numerically. The red dashed lines are guide 
to the eye for the positions of the Weyl nodes. The parameters in (I)-(VI) are in one to one correspondence to the six cuts 
labeled by (I)-(VI) in Fig.S |T] Note that the ir-axis in Fig.S |T] is m z /mo instead of m z . 


In the SCBA, which specializes to weak disorder, the self-energy D(Ep) can be obtained by solving the self-consistent 
integral equation, 


W 2 a 2 (' 

S (^) = nr^/ d 3 k[(E F +i0+)I-H o -E]-\ (6) 

967r d J kGBZ 

where the k integral is over the first Brillouin zone (BZ) and the lattice constant a = 1. The solution of the self-energy 
can then be used to obtain the renormalized band dispersion in the SCBA. In solving the integral equation (|5|l . we 
set t x ,y,z = l,wo = 2.1 and Ep = 0 as in the main text. It is straightforward to verify that the Fermi level remains 
unrenormalized, i.e. Ep = 0. 
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In Fig.S [2l we plot the calculated SCBA band dispersion for different disorder strength W = {0, 2,3,4} and in 
systems with different bare mass m z = { — 1.2, —1, —0.8,0.8,1,1.2} sitting near the QAH-WSM (fh z /mo = — t~/m,o) 
and the WSM-NI ( fh z /mo = t z /mo) phase boundaries as indicated by the corresponding roman numerals in Fig.S [T] 
Figs.S[2}II-III) show that for —fh z < t~, the Weyl nodes located at kP z = ± arccos(m z /t z ) are close to the BZ boundary 
and will move to and annihilate at the BZ boundary (fcg = 7r) with increasing W. As a result, the WSM enters the 
QAH phase by opening a nontrival (inverted) band gap with fh z < —t z that continues to increase with disorder [see 
Figs.S EDI)] - In contrast, for fh z < t z , Figs.S [2{IV-V) show that the Weyl nodes repel each other near the BZ center 
(k z = 0) as the disorder effectively increases their separation by increasing the renormalized mass fh z . Interestingly, 
in Fig.S [2]( VI) where fh z > t z , a pair of Weyl nodes is nucleated at zone center with increasing disorder strength, 
which gives raise to a normal band insulator to WSM transition. 

Furthermore, we can obtain a closed form of £ and gain more analytical insights by neglecting £ on the right side 
of EqO which amounts to the Born approximation. Expanding Hq to quadratic order in k , we obtain 


m z 


m z 


W 2 a 2 1 
967T 2 mo 



m 0 

- t z kl - m z 



dk z . 


This result shows that for the parameter regime studied in the main text, the disorder-induced correction to fh z always 
has the same sign as —mo. The band gap therefore becomes more inverted with increasing disorder strength W. In 
particular, for a WSM with — t z < m z < t z , the pair of Weyl nodes move to the BZ boundary where fh z approaches 
— t z and annihilate. The WSM will then open up an inverted topogolical gap (with rh z < —t z ) that further increases 
with the disorder strength W. If the clean system is a normal band insulator with m z > t z , the normal band gap will 
close (when m z = t z ) with increasing W by nucleating a pair of Weyl nodes at the zone center that repel each other 
upon further increasing W. 
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